Dynamics of Dissipative Quantum Systems — from Path Integrals to Master Equations* 

Joachim Ankerholcjj] 

Physikalisches Institut, Universitat Freiburg, Hermann-Herder-Strafle 3, 79104 Freiburg, Germany 

The path integral approach offers not only an exact expression for the nonequilibrium dynamics 
of dissipative quantum systems, but is also a convenient starting point for perturbative treatments. 
| An alternative way to explore the influence of friction in the quantum realm is based on master 

, equations which require, however, in one or the other aspect approximations. Here it is discussed 

under which conditions and limitations Markovian master equations can be derived from exact path 
integrals thus providing a firm basis for their applicability. 
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I. INTRODUCTION 

Quantum systems coupled to a heat bath environment can be found almost everywhere in physics and chemistry. 
Recently, challenging problems e.g. in realizing quantum computers, in manipulating Bose-Einstein condensates, or 
in understanding low temperature dynamics in molecular systems have triggered a substantial amount of research. 

The description of quantum dissipative processes reaches back to the late 1960s 0, when it was mainly concerned 
with weakly damped quantum optical systems. The conventional framework relied on quantum Langevin or quantum 
master equations. A breakthrough beyond the limitations of the weak coupling approach was made in the early 1980s 
, t | 0- Path integral techniques were shown to be powerful means to formulate the reduced dynamics of dissipative 
C$ . quantum systems for all damping strengths, temperatures and bath memory times. While applications, particular 
in condensed phase systems, have proven the advantage of this approach, e.g. the non-exponential decay of low 
temperature correlation functions, it has rarely been used to put the known master equations on a firm basis and to 
derive new master equations for strong damping/low temperatures. In fact, it turned out that the quantum stochastic 
process is strongly non-Markovian and intimately depends on the initial correlations between system and bath 0. 
As a consequence, "simple" master equations do in general not exist and the known results can be derived only in 
an approximate way What has been essentially unexplored so far, namely the range of strong friction, has been 
analyzed recently in detail 0, . It has been shown that in this limit quantum noise is squeezed and that the Wigner 
transform of the reduced density obeys a quantum Fokker Planck equation Q. In the low temperature quantum 
domain the latter one can be reduced to a generalization of the classical Smoluchowski equation 6] . Since the path 
integral formulation is exact in the system-bath coupling and includes also correlated initial states, it is basically the 
CNJ , only way to derive these findings consistently. 

0^ ' Thus, we are now in a position to give a more or less comprehensive account about the validity and limitations of 
Markovian master equations. The following study is a brief attempt in this direction. It is not intended, however, to 
include all developments and to consider the most general case. Instead, to keep things as transparent as possible I 
restrict myself to (i) a single one-dimensional continuous system with a well defined ground state coupled to a single 
heat bath where (ii) the bath cut-off frequency is assumed to be the largest frequency scale. Further, I look only for 
Markovian master equations with (iii) time evolution generators independent of time and initial preparation. It turns 
out that already within this frame the analysis is rather subtle and reveals most of the complexity we encounter in 
'■^J ' describing dissipative dynamics. 

In Sections ^ and I1VI I briefly introduce the path integral approach for dissipative quantum systems. Based on 
this in Sec. II VI some general conditions for the existence of Markovian master equations are derived which are then 
specified in Sec. for weak and in Sec. IVII for strong damping, respectively. 



o 



- 1—1 

X 



II. PATH INTEGRALS FOR DISSIPATIVE QUANTUM SYSTEMS 

The description of classical Brownian motion in terms of generalized Langevin equations or, equivalcntly, in terms of 
Fokker-Planck equations for phase space distributions has a long tradition 8]. In contrast, the inclusion of dissipation 
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within quantum mechanics in a non-perturbative way has been established only since the early 80s |2( ■ In the standard 
formulation one starts with a system+heat bath formulation 

H = H S +H K + Hi (1) 

where the total Hamiltonian contains a system, a reservoir (heat bath), and a system-bath interaction part, respec- 
tively. The dynamics of the total density matrix starting at t = from a general initial state W(0) evolves according 
to 



W(t) = exp(-iHt/K) W(0) exp(iHt/h). (2) 

The crucial point is now, that dissipation is not a priori inherent in the system, but arises only if one looks on the 
effective impact of the bath degrees of freedom within a reduced picture g(t) = ti^{W(t)}. The Gaussian statistics 
of the heat bath is modeled by a quasi-continuum of harmonic oscillators bilinearly coupled with the relevant system 
degree of freedom. Although the interaction between each bath degree of freedom and the system is supposed to be 
weak, the overall impact of the reservoir may cause also strong friction. Along this reasoning two steps need to be 
done: an appropriate formulation to arrive at a reduced description has to be found and the initial state has to be 
defined. 

The only non-perturbative way to deal with the elimination of the bath degrees of freedom is to apply the path 
integral approach. Denoting the degrees of freedom in coordinate space by x for the bath and q for the system, the 
coordinate representation of (J2J) follows as 

(q f ,Xf\W(t)\q' fl x' f ) = J dq i dq' i dxidx!G(qf,Xf,t,qi,£i) 

{q u 2i|W(Q)|<«) G(q' f ,x' f ,t, q\, 31)* (3) 
where the * means complex conjugation. The transition amplitudes on the rhs are expressed as path integrals, e.g., 

G(q f ,x f , t, q^xi) = J V[q] V[x\ e iS ^' h (4) 

with the total action S = Sg + Sr. + Si according to the three parts of the Hamiltonian (J2J. The sum goes over 
all paths running in time t from g(0) = qi,x(0) = x% to q(t) = qf,x(t) = Xf. Switching to a reduced description is 
achieved by taking the trace over the bath degrees of freedom 

Q(qf,q' f ,t)= J dx f (q f ,x f \W(t)\q' f ,x f ). (5) 

To carry out all integrations over the bath degrees of freedom in (0J explicitly, the initial state must be specified. 



III. INITIAL STATE AND INFLUENCE FUNCTIONAL 



In the ordinary Feynman Vernon theory [2| the initial state is assumed to be a factorizing state W(0) = 
Ps(0) Z^ 1 exp(— /3-Hr) (^r is the partition function of the bath) so that each one, system and equilibrated bath, 
lives in splendid isolation at t = 0. While this assumption may be justified in the weak friction/high temperature 
range, it definitely fails for moderate to strong dissipation or lower temperature. Even the Langevin equation is not 
regained in the classical limit, but differs by initial boundary terms that may persist up to long times. A realistic 
initial state reflecting the experimental situation is thus a correlated one described by Q 

vno) = vE °y~ PH °i ( 6 ) 

i 

where Zp is the partition function of the total system and the operators O l s , Og act onto the system degree of freedom 
only and prepare a nonequilibrium state. In the sequel I focus on the case where the preparation operators depend 
exclusively on coordinate and refer to 3] for the generalization. As an example think about a position measurement 
with a Gaussian slit, in which case the preparation operators are Gaussian weighted projection operators onto position. 
Representing the statistical operator in © as path integrals in imaginary time we have 



( qi ,Xi\W(0U, Si) = Z/ A(g 4s ^) J V[q]V[i\ e~ s ^l n 



(7) 
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with the preparation function A(-) being the coordinate representation of the preparation operators in JHJ and the 
total Euclidian action S = Sg + <Sr + Si. Paths contributing to J7J) run in imaginary time in the interval h(3 from 
q(0) = q i: 1(0) = Xi to q{h{3) = q[, x(k/3) = x[. 

The integrations over the bath degrees of freedom in © can now together with Q be performed exactly due to 
the harmonic nature of the bath. After some tedious algebra one ends up with the coordinate representation of the 
reduced density matrix 



Q{qf,q'f,t) = J d^dq- J(qf,qf f , \(qi,(£). (8) 
The propagating function J(-) is a threefold path integral over the system degree of freedom only 

J(q f ,q' f , t, q liq [) = Z- 1 f V[q] V[q'] V[q] j^M Ml* (9) 



with Z = Zp/Zn. The two real time paths q(s) and q'(s) connect in time t the initial points qt and q[ with the fixed 
end points qf and q'p while the imaginary time path q(a) runs from q^ to q[ in the interval fi(5. The contribution 

of each path is weighted with an effective action T,[q, q',q] — Ss[q] — Ss[q'] + iSs[q] + i<P[q, q' , It consists of the 
actions of the bare system in real and imaginary time, respectively, and additional interaction contributions (influence 
functional) 

4>[q,q\q] — ~* / dr / daq(r)K(— ir + ia) q(a) +i / dr— q(r) 2 
Jo Jo Jo * 



rh/3 rt 

+ I dr / ds K*(s — ir) q(r) x(s) ~ Mrij(s) x(s) 
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t 

ds x(s) 



M I du / -f(s — u)r(u) / du K'(s — u) x(u) 
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(10) 



Here, for convenience, we have introduced sum and difference real time paths, namely, r(s) = [q(s) + q'(s)]/2 and 
x(s) — [q(s) — q'(s)]. Further, the effective impact of the bath is controlled by the damping kernel 

J 7T smh(ujh(j/2) 

where 9 = s — ir, < s < t, < r < hf3 and I(u>) is the spectral density of the heat bath. For an interaction term 
between system and bath of the form H\ = q£n(6) with £r(0) = cx (coupling constants c) one shows that the damping 
kernel is basically the bath autocorrelation function, i.e. K (9) = (£r(0) Cr(0))r/?i. In real time K(s) — K'(s) + iK"(s) 
is related to the macroscopic damping kernel 

2 f°° du; I(u) 

l(s) = — / — -^cos(w«) 12 

via K"{s) — ^-d , y(s)/ds and K'(s) — > Mj(s)/frf3 in the classical limit loJi(5 — > (M is the mass of the Brownian 
particle). The term with = lim^^o h(3K(0) in the first line in i|10[l gives a potential renormalization in the 
Euclidian action due to shifts of the minima of the bath oscillators by the coupling to the system. The corresponding 
renormalization in the real-time part of the action has already been incorporated in the form it is written in i|10fl • 

Apparently, (p[q, q',q] contains the bath induced non-local interactions between the various system paths. Particu- 
larly, it turns out that the last term (with K') acting as a Gaussian weighting factor for the quantum fluctuations x(s) 
suppresses coherences and gives rise to relaxation, while the second last term depending on K"(s) leads to dephasing 
during the dynamics. The influence of initial correlations between system and bath is contained in the term with K* 
coupling real and imaginary time motion. 

While in JHJl the imaginary time paths describe the initial state, the two real time paths govern the dynamics 
of the reduced system. Accordingly, the distribution of end-points of the former and starting points of the latter 
qi,q[ are weighted in (JHJ also by the preparation function A(-). In the limit t — > one has J(qf,q'j : ,t 1 qi,q' i ) — > 
Qpiqi.q'i) 5{q f - qi)8{q' } - q[) so that 

Q{Qf,rff, °) = Qeiquq'i) Kiirfi) ( 13 ) 

with the reduced equilibrium density matrix Qp(q, q') — Z^ 1 (g|trR exp(— (3H)\q'). In fact, this formulation reproduces 
in the classical limit the generalized Langevin equation. 
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The nonequilibrium time evolution of a dissipative quantum system is governed by (JSJ together with l|13(l . The 
good news is that this path integral expression is exact, also in the system-bath coupling. The bad news is that 
its evaluation is even numerically feasible only for specific systems and certain ranges in parameter space. Since 
the propagating function is highly oscillatory, numerical algorithms inevitably become unstable for sufficiently long 
times. In this situation "simple" time evolution equations are of great practical importance and the path integral 
formulation provides a very convenient basis for deriving them. The severe problem, however, is the non-locality, in 
time, of the influence functional; the reduced time evolution in the time interval s£ [t',t] C [0,t] is affected by the 
history of the dynamics for < s < t' and particularly by the initial correlations between system and bath at s = 0. 
Even worse, it has been shown that different initial densities W(0) of the full compound may lead to identical reduced 
initial densities p(0), but give rise to quite different quantum stochastic processes due to different initial correlations 
3j. Hence, in general a time evolution equation g = C g of the reduced density with a generator C independent of 
time and independent of the initial preparation does not exist. 

I mention here that an alternative procedure to eliminate the bath degrees of freedom has been developed by 
invoking projection operator techniques 0,0. This way generalized equations of motion in form of integro-differential 
equations for the reduced density have been derived. It turns out, however, that this approach is convenient only in the 
weak damping/high temperature range. In contrast, the path integral formulation offers an elegant starting point for 
perturbative treatments, e.g. semiclassical approximations, perturbation theory a la Feynman etc., for all temperatures 
and spectral bath densities. Moreover, its close relation to the formulation of classical statistical mechanics has allowed 
to adopt and extend numerical techniques like e.g. Monte Carlo simulations to quantum dissipative systems. 



IV. GENERAL CONDITIONS FOR MARKOVIAN MARKOVIAN MASTER EQUATIONS 



As the time evolution of the reduced density matrix cannot be cast into the form g = C g exactly, one may wonder 
under which conditions at least approximate time evolution operators £, independent of time and the initial state, 
can be derived from the expression JSJ. Before I turn to specific cases, in this section, we first want to formulate some 
general conditions. Accordingly, we analyze the influence functional (|10|l with respect to its non-Markovian nature 
and with respect to the friction induced entanglement of the bare time evolution operators. 

For this purpose let us consider the real-time part of the damping kernel in the limit of uj c — > oo. One has 



K(t) 



1 



(hp) 2 smh 2 {nt/h[3) 



+ iMjS(t). 



(14) 



Hence, while K"(t) is local in time, K'(t) is not. In fact, for h(3 — > oo it only decays algebraically K'(t) cx 7/i 2 
so that strong non-local quantum fluctuations in the influence functional l|10(l become important and lead e.g. to 
non-exponential long time tails for zero-temperature correlation functions Q|. 

One concludes that a time evolution operator, local in time, can only exist on a coarse grained time scale s 3> 
hf3,l/uj c . In fact, a coarse graining procedure is only meaningful if the time scale on which the relevant dynamics 
takes place, e.g. the relaxation time scale t r , obeys 



t r > hp, l/w c 



(15) 



Compared to the typical time scale of the bare system time evolution 1/u)q with ujq e.g. its ground state frequency, t r 
can be very large. For T — > the thermal time scale exceeds all other time scales and the above condition can never 
be fulfilled. 

The immediate consequence of <|15fl is that the influence functional becomes local on the coarse grained time 
scale. This is fine, but it is not all we need. Namely, the bath mediated coupling between forward and backward 
time evolution and also between the time evolution and the initial state — accounted for by the x dependence of the 
influence functional (|10f> — is a genuine quantum effect, so that the x paths can be interpreted to describe quantum 
noise. For anharmonic potential fields V(r — x/2) the effective action contains also anharmonic x dependent terms 
meaning that quantum noise is in general non-Gaussian. In the spirit of a classical Kramers-Moyal expansion ||| one 
would expect that a time evolution generator L of the reduced dynamics p(x/,rf, t) must be represented as an infinite 
power series in Xf and d/drf (leading in the Wigner transform to derivatives d/dp, d/dq) with an infinite number of 
diffusion coefficients — a hopeless situation. For further exploring this point it is useful to introduce the dimensionless 
quantity 



1 



ds du K' (s — u) 



2t r 7 
UlQ hf3 



(16) 
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This parameter is obtained by considering the last term of l|10(l . which acts as a Gaussian weight for the x-paths, for 
times of the order oft r , putting x(s) — const., exploiting l|14|). and performing the coarse graining. Accordingly, the 
typical size of \x\/ y/h/Mcuo is approximated to be of the order 1/t/k and quantum noise tends to be small in domains 
of parameter space where k» 1. 

In the overdamped range j/ujq ^ 1 relaxation takes place on the scale t r = j/ll>q so that n — (7/^0) j/uoh/3 is 
large under l|15|l . For sufficiently strong damping the Gaussian weighting factor in the influence functional then causes 
quantum noise to be squeezed by friction. This allows for a semiclassical type of approximation with generators C 
containing at most second order diffusion coefficients (ranges I and II in Fig. In the underdamped case j/u>q < 1 
the dissipative system equilibrates on the time scale t r ~ I/7 so that l|15l) reads 7?i/3 <C 1 and n — 2/u>ofif3 is 
large as well. We this way recover the well-known fact that in the high temperature range luqHP <C 1 quantum 
fluctuations are small and quantum nonequilibrium dynamics happens to be close to classical nonequilibrium dynamics 
(range III in Fig. The coarse graining condition in the underdamped regime 7?i/3 <C 1 covers also the range IV 
in Fig. ^ There, quantum noise is non-Gaussian and not small at all but the entanglement due to friction can 
be treated perturbatively. To see this we estima te the contribution of kinetic terms in the bare real-time actions 
to be of order M{x/Ti(3) 2 (t r /h) — (x/ y/fi/Muj^) 2 l/[^ujo(hp) 2 ], while the x 2 term in the influence functional is 
(xj y^h/Mujo) 2 k — (x I ^/fif Mujq) 2 2/ ujoh(3 . Thus in the domain where the condition l)15|l is obeyed for j/ojo < 1, i.e. 
"fhp <C 1, friction terms are supposed to be much smaller than bare kinetic contributions. 

To summarize this discussion we can expect to derive from the exact path integral expression JHJl approximate time 
evolution operators £ on a coarse grained time scale only in the domains I-IV of Fig.^ In the sectors I— III quantum 
noise is essentially Gaussian — either due to strong friction (I) or due to high temperature (II, III) — while in sector IV 
it is non-Gaussian, but the coupling it induces between the foward and backward time evolution of the bare system is 
weak. In the remaining parts of parameter space there is no reduction in a simple way possible and one has to work 
either with the full path integrals or master equations with complicated non-local integro-differential operators. In 
the sequel I will focus on the ranges III, IV for weak damping and I, II for strong friction in detail. 



V. WEAK DAMPING REGIME 



It is well-known that for sufficiently high temperatures and/or sufficiently weak friction time evolution equations for 
the reduced dynamics can be derived [lj. This has been done already in the late 60s motivated by the experimental 
progress for quantum optical devices. Here, I discuss how these results in the underdamped range can be regained 
from the exact path integral expression JSJ). 

In case where fluctuations tend to be small uioh(3 <C 1 (range III in Fig. ^) master equations have been explic- 
itly derived for harmonic systems. From (JSJ an equation of motion for g(xf,rf,t) is obtained which, in operator 
representation, coincides with the well-known Agarwal equation |lOl |: g(t) = C A g o{t) with 



c Ag e = ~[h , e ] - l ^[ q , {p, e }] - (r • , [<i- [<i- <■]] ( i < ) 



where H — p 2 /2M + Muj^q 2 /2 and [•, •] ({•, •}) denotes the commutator (anticommutator) . The equilibrium position 
variance is approximated to read 

(^-T^coth^V (18) 



Mw 

Since for a harmonic system quantum noise is always Gaussian, the above finding can be generalized to the wider 
range jhp <§; 1 (ranges III and IV in Fig.QJ. The result is an additional p-q diffusion term, i.e. 

C cxt q = C Ag Q + [p, [q, £>]] ■ (19) 

The form of this additional diffusion coefficient for the harmonic case has been first derived from the path integral 
approach in Q 



' I (Mufa 2 )?-^) (->()! 



with (p 2 )f3 the equilibrium momentum variance of a harmonic oscillator. For the parameter range considered here 
7^/3 <C 1 (u)q <C cj c ) one has 



b?)p w M 2 l4 (q 2 ) + ^ [ V , (i + Wc ^/27r) + C] (21) 
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FIG. 1: Sectors in parameter space where for u> c 2S> Wo,7 various types of master equations can be derived I-IV. The thin 
horizontal line separates the overdamped from the underdamped range, the thick solid line defines the range S> 1 in the 
overdamped region, the dashed line defines yhp <C 1, and the dotted vertical line specifies uioh/3 < 1 in the underdamped 
sector. See text for details. 
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where ip(x) is the psi function and C w 0.577 Euler's constant. Accordingly, we infer from (|20|l in ranges III/IV 

D pq « -2 [-0 (1 + w c ^/3/27r) + C] . (22) 



7T 



This Z? pg tends to its classical limit D P9 = for ui c %(3 — > with fixed w c 3> t^o- From the Wigner transform of £ ox t 
we then regain the classical Fokker-Planck operator [8j. For finite but still small u}qU(3 <C 1 the influence of D pq is 
supposed to be small only, if additionally u> c h(3 <C 1, which gives us the precise validity of the Agarwal equation l|17(l 
in range III. For lower temperatures, sectors III (partially) and IV, the p-g diffusion becomes important and £ ext with 
D pq as in 1)22(1 defines a master equation derived by Haake and Reibold [Ij. Within the Nakajima-Zwanzig projection 
operator technique 9] an identical equation has been obtained within Rcdficld theory [12}. In the limit ui c hj3 3> 1 the 
expression (|2~2"|l simplifies to leading order to D pq sa (7/71- ) \og{uj c Ti (3 /2-k). 

So far I have discussed only master equations for harmonic systems. In case of anharmonic systems and within 
the range ujqTiP *C 1 in principle a semiclassical type of approximation applies so that the potential field V(r — x/2) 
can be expanded up to second order in x. This way, as mentioned above, it is shown that the x— paths describe 
effectively white Gaussian noise in the reduced system. Diffusion related terms thus appear in C at most in second 
order and the Wigner transform of C into phase space contains at most second order derivatives in the phase space 
variables leading to a generalized Fokker-Planck operator. However, according to the previous section, in the wider 
range III/IV a weak damping perturbation theory of the bath induced entanglement is feasible anyway. Exploiting 
this very fact, the noncquilibrium quantum dynamics can be cast, after an additional coarse graining s 3> 1/wo, into 
the form of the famous Lindblad equation 13] 

m = ~[h\ g (t)} + ^ E m*)> L \\ + L <> emh ■ ( 23 ) 

1 



Here, H' is a hermitian operator which must not necessarily coincide with the bare system Hamiltonian, and the 
Li describe the effective influence of the specific heat bath. In case of a harmonic system and u sing creation and 
annihilation operators the resulting master equation is also known as the Haake- Weidlich equation |l4| . 

Before I turn to the strong friction range two remarks are in order here. First, all time evolution operators C 
specified in (|17fl and l|19l) are not of Lindblad form (|23|l and thus do not respect complete positivity. This is not 
astonishing since they are derived based upon a coarse graining procedure. Their validity is restricted to a certain 
subspace where transient components have died out, while Lindblad theory itself requires a Markovian master equation 
for all times. The path integral approach reveals explicitly that from a physical point of view this requirement is 
not necessary (see also |l5| ) and even does not reflect the exact nonequilibrium dynamics. Second, master equations 
have also been derived starting with factorizing initial states (Feynman- Vernon theory), see Sec.[H] e.g. the Caldeira- 
Leggett master equation |l6j . The problem there is that the factorizing initial state itself makes only sense in the 
weak damping limit meaning that a weak friction perturbation theory within Feynman Vernon theory is somewhat 
inconsistent. Particularly even the known Fokker-Planck operator cannot be re-derived in the classical limit |l7j . 
Hence, the Caldeira-Leggett or related master equations cannot be obtained from @ Q|. 



VI. STRONG FRICTION RANGE 



In classical physics the strong friction domain is well-known as the Smoluchowski limit. Its characteristic property 
is that one has a separation of time scales between fast equilibration of momentum and slow relaxation in position. 
This way the Fokker-Planck equation for the phase space distribution can be reduced to a Smoluchowski equation 
for the marginal distribution in position space. For quantum dissipative systems one would think that strong friction 
makes the reduced system to behave effectively more classically so that the complicated path integral expression © 
simplifies considerably. That this is indeed the case has been analyzed in detail only recently 0, Q and I have given 
the general argument already in Sec. II VI In particular, strong dissipation has quite a different influence on position 
and momentum. 

To see that in detail we define a typical damping strength as 

7 ee 7 (0)= lim M (24) 

with 7(0;) the Laplace transform of the classical damping kernel ^(t) l|12fl . For instance, in the ohmic case I(u>) = Mjlu 
and also for the more realistic Drude model I{lu) = Mjluuj^/(u! 2 + ui^) one finds 7 = 7. Given a typical frequency loq 
of the bare system, e.g. the ground state frequency, by strong damping we then mean 

J_» ft/3, 1,1. (25) 

Wn UJ C 7 
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In other words, we assume the time scale separation well-known from the classical overdamped regime |8j and extend 
it to the quantum range by incorporating the time scale for quantum fluctuations %0. Correspondingly, according to 
the discussion in Sec. II VI we consider the dynamics (jBJ on the coarse grained time scale s hp, — and a ^> 
The consequences are threefold: (a) the strong friction suppresses non-diagonal elements of the reduced density matrix 
during the time evolution; (b) the real-time part K(s) of the damping kernel becomes local on the coarse grained time; 
(c) initial correlations as described by the first term in (|10ll survive for times of the order 7/<x>q so that factorizing 
initial states cannot be used. 

Following the above simplifications the path integral formulation now allows for a perturbative treatment in the 
strong damping limit. According to Sec. the strategy is then to evaluate the path integrals in the sense of a 
semiclassical approximation by assuming self-consistently that non-diagonal elements, i.e. x — q(hf3) — q(0) and x(s) 
dependent terms, remain small during the time evolution. Hence the effective action E[r, x,q] is expanded up to 
second order in % of the imaginary time paths and in the excursions x(s) of the real-time path integrals. Doing so we 
take sufficiently smooth potentials for granted. It is worthwhile to note that we do not need to restrict the value of 
jhp meaning that the strong damping analysis covers a broad temperature range from the classical (jhp <C 1 , range 
II in Fig.^) to the deep quantum domain ("fhf3 3> 1, range I). 

Before I turn to the dynamical case, it is instructive to consider the thermal equilibrium only. The corresponding 
imaginary time path integral can approximately be solved for an anharmonic potential and the result is 

gp (x, f) = - e-P V ^)-(P 2 ) * 2 /^ 2 e A/3[/3V'( F ) 2 /2-3V"(r)/2] (2g) 

where f = [q{%0) + q(0)]/2 and Z = J dqg p (0, q). Further, 

Apparently, A measures the typical strength of quantum fluctuations in position space. In case of Drude damping with 
a high frequency cut-off ui c both A and (p 2 ) can be expressed in terms of ^> functions. Then, for high temperatures 
jh/3 « 1 we find A ps Ti 2 j3/12M and (p 2 ) ~ AI//3. The friction dependence appears as a genuine quantum effect 
for lower temperatures and for jhf] ^> 1 one has A ss (h/M^Tr)\og('jfi(3/2ir) and (p 2 ) « (Mh'y/ir) log(w c /7). With 
increasing 7 the strong squeezing of quantum fluctuations in position induces enhanced quantum fluctuations in 
momentum, thus suppressing non-diagonal elements in the density matrix. Interestingly, the probability distribution 
is Gaussian in x, i.e. its Wigner transform (x/h — > p) Gaussian in momentum, even at low temperatures. Anharmonic 
corrections in x to the exponent are at most of order I/7 2 . Hence, for strong friction the equilibrium density consists 
of a part which in phase space takes the form of a classical distribution, however, with an h dependent (p 2 ) and a 
part with A dependent quantum corrections. 

In a similar way the dynamics can be treated. In essence, since deviations from diagonality x(s) = q(s) — q'(s) remain 
small during the time evolution, they run effectively at each instantaneous mean position r(s) = [q(s) + q'(s)]/2 in a 
harmonic force field V"(r)x. Exploiting also the sluggish motion of r (s), this allows for an analytical solution which 
eventually leads to a time evolution equation for g{qf, q'f,t). After switching to classical phase space {xf, ry} — > {p, g}, 
i.e. g(xf,rf,t) — » W(p,q,t), one arrives at 



5»w> - {5 Pito+'H'! -££+■»<■•>>!, 

+3^ [1/0 + A V"(q) - {p 2 )/M] \ W(p, q,t). (28) 

Here, we have introduced an effective potential V e s = V + (A/2) V". The first line on the r.h.s coincides with a 
classical Fokker-Planck operator in an effective force field, the second line describes quantum mechanical coupled 
p-q diffusion, a necessary ingredient as already seen for the weak damping case in the previous section. In the high 
temperature limit 7?i/3 — > the quantum Fokker-Planck equation (QFP) tends to the classical Kramers equation y| • 
Small but finite jhp <C 1 means in the overdamped region also uofiP <C 1 so that in case of a harmonic potential 
the QFP becomes identical to the master equation gained by Haake and Reibold specified through (|I9|I and l|2U|l . 
However, while this known master equation is restricted (for 7/c^o > 1) to the range 7?i/3 <C 1, the above QFP is valid 
for all and for anharmonic potentials as well. It is also not of Lindblad form due to the coarse graining procedure 
on which its derivation is based. Of course, the equilibrium solution to l|28|l is given by the Wigner transform of 
Let us briefly touch the question about higher order diffusion terms to l|28|l . For harmonic systems they do not occur 
so that the QFP is in this sense exact [Tg ■ In case of anharmonic potentials they result from non-Gaussian quantum 
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TABLE I: Collection of approximate master equations that have been derived within the exact path integral approach in certain 
domains of parameter space. 



Master equation 


Validity 


harmonic 


anharmonic 


Remarks 


Agarwal 


in 


yes 


no 


uj c h(3 < 1 


Haake-Reibold 


II, III, IV 


yes 


no 


derived also for uj c > uia 


Haake-Weidlich 


III, IV 


yes 


no 


coarse graining t 3> 1/loq 


Lindblad 


III, IV 


yes 


yes 


coarse graining t ^> 1/uo 


QFP 


I, II 


yes 


yes 


can be reduced to QSE 



a Range of validity according to Fig. 



fluctuations attributed to higher than second order derivatives in V(q). A rough estimate shows that anharmonic 
terms in xj (leading to higher than second order derivatives in p) in the crucial low temperature range "fhf3 3> 1 are 
of order l/[h'-/ 3 ^ 2 \og(ui c /j)] compared to the leading terms. 

With the QFP at hand, we are now able to follow the reasoning of classical mechanics and reduce it to an equation 
in position space- the quantum Smoluchowski equation (QSE). For this purpose we employ the projection operator 
techniques invoked in Along these lines one introduces the operators V = fpip) J dp and Q = 1 — V where fp{p) 
is the normalized momentum distribution in equilibrium according to (|26[) . The next steps are straightforward and 
not presented here in detail. After some algebra one arrives to order A/7 3 at an equation for the position distribution 
n (q, t) = f- 1 VW of the form 

§l n ^ t] = W I { 1 + Wf + Av "' ] } L ^ n ^ <} • (29) 

Here, to leading order the time evolution is determined by 

i Q SE = Kff + d/dq[l/P + AV"] (30) 

derived first in 0. Inertia corrections appear in curly brackets: A classical (A = 0) correction oc V" shows up, 
while quantum fluctuations enter through third and forth order derivatives of the potential. Overdamped quantum 
Brownian motion in position space thus becomes much more sensitive to the details of the potential profile at lower 
temperatures. 

VII. CONCLUSIONS 

Based on the exact path integral approach I have explored to what extent the nonequilibrium dynamics of a 
dissipative quantum system can be described by a time evoution operator independent of time and initial preparation. 
The main results of this study are summarized in Table 1 together with Fig. ^ The Agarwal equation, the Haake- 
Reibold, and the Haake-Weidlich master equations are restricted to harmonic systems, where only the Haake-Reibold 
result is also applicable for strong damping/high temperature. The Lindblad equation can be used in the weak 
damping domain for general potential fields but includes an additional coarse graining on the time scale of the bare 
system time evolution. For strong damping and high as well as low temperatures only the recently found quantum 
Fokker-Planck/Smoluchowski equation approximates the exact dynamics. Combining both, Lindblad's equation and 
the latter one, allows to conveniently treat dissipative quantum systems in a considerable range of parameter space. 
To obtain corresponding generators £, a coarse graining procedure is required which renders either a semiclassical 
type of approximation due to squeezed Gaussian quantum noise or a perturbative treatment due to weak bath induced 
entanglement to be successful. 

The area in Fig. ^ where no reduction to a simple time evolution equation is possible is a real challenge for 
further developments. One direction in which progress has been made uses Hubbard-Stratonovich transformations to 
disentangle forward and backward paths in the influence functional (|10|) at the expense of new auxiliary fields. Along 
these lines the reduced dynamics has been expressed in terms of stochastic, nonlinear Schrodinger equation with 
complex noise j2Cj. In the same spirit, the quantum stochastic process has been described to involve discontinuous 
jumps 23]. So far, however, the applicability of related methods is restricted to specific systems and/or certain regions 
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in parameter space. Efficient algorithms to solve the exact path integrals JHJ have been successfully developed and 
applied, e.g. real-time Monte Carlo schemes 21] or the so-called QUAPI-method [2^, but the broad range of long 
times, intermediate damping, low temperatures, and long bath memory times is still out of reach. Strong efforts to 
overcome this lack are being on the way. 

I am indebted to H. Grabert and P. Pechukas for many fruitful discussions and suggestions. 



[1] For a review see: F. Haake in Quantum Statistics in Optics and Solid-State Physics, G. Hohler (ed.), Springer Tracts in 

Modern Physics 66 (Springer, Berlin, 1973). 
[2] U. Weiss, Quantum Dissipative Systems, (Singapore, World Scientific, 1999). 
[3] H. Grabert, P. Schramm, G.-L. Ingold, Phys. Rep. 168, 115 (1988). 
[4] R. Karrlein and H. Grabert, Phys. Rev. E 55, 153 (1997). 
[5] M. Thorwart, M. Grifoni, and P. Hanggi, Phys. Rev. Lett. 85, 860 (2000). 
[6] J. Ankerhold, P. Pechukas, and H. Grabert, Phys. Rev. Lett. 87, 086802 (2001). 
[7] J. Ankerhold, Euro. Phys. Lett. 61, 301 (2003). 
[8] H. Risken, The Fokker Planck Equation, (Springer, Berlin, 1984). 

[9] S. Nakajima, Progr. Theor. Phys. 20, 948 (1958); R. Zwanzig, J. Chem. Phys. 33, 1338 (1960). 
[10] G. S. Agarwal, Phys. Rev. 178, 2025 (1969). 
[11] F. Haake and R. Reibold, Phys. Rev. A 32, 2462 (1985). 

[12] A.G Redfield, Adv. Magn. Reson. 1, 1 (1965); K. Blum, Density Matrix Theory and Applications, (Plenum, NY, 1981). 
[13] G. Lindblad, Commun. Math. Phys. 48, 119 (1976). 
[14] W. Weidlich and F. Haake, Z. Phys. 185, 30 (1965). 
[15] P. Pechukas, Phys. Rev. Lett. 73, 1060 (1994). 

[16] A.O. Caldeira and A.J. Leggett, Physica A 121, 587 (1983); B.L. Hu, J. P. Paz, and Y. Zhang, Phys. Rev. D 45, 2843 
(1992). 

[17] P. Hanggi, Z. Phys. B 31, 407 (1978). 

[18] P. Talkner, Ann. Phys. (NY) 167, 390 (1986). 

[19] J. L. Skinner and P.G Wolynes, Physica A 96, 561 (1979). 

[20] L. Diosi, N. Gisin, and W.T. Strunz, Phys. Rev. A 58, 1699 (1998); J.T. Stockburger and H. Grabert, Phys. Rev. Lett. 
88, 170407 (2001). 

[21] R. Egger and C.H. Mak, Phys. Rev. B 50, 15210 (1994); C.H. Mak, R. Egger, and H. Weber-Gottschick, Phys. Rev. Lett. 
45, 33 (1998). 

[22] D.E. Makarov and N. Makri, Chem. Phys. Lett. 221, 482 (1994). 
[23] H.-P. Breuer and F. Petruccione, Phys. Rev. Lett. 74, 3788 (1995). 



